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ABSTRACT 

We introduce a novel method to measure the masses of galaxy clusters at high redshift selected from 
optical and IR Spitzer data via the red-sequence technique. Lyman-break galaxies are used as a well 
understood, high-redshift background sample allowing mass measurements of lenses at unprecedented 
high redshifts using weak lensing magnification. By stacking a significant number of clusters at 
different redshifts with average masses of ~ 1 — 3 x lO^^M©, as estimated from their richness, we can 
calibrate the normalisation of the mass-richness relation. With the current data set (area: 6 deg^) 
we detect a magnification signal at the > S-cr level. There is good agreement between the masses 
estimated from the richness of the clusters and the average masses estimated from magnification, 
albeit with large uncertainties. We perform tests that suggest the absence of strong systematic effects 
and support the robustness of the measurement. This method - when applied to larger data sets in the 
future - will yield an accurate calibration of the mass-observable relations at z > 1 which will represent 
an invaluable input for cosmological studies using the galaxy cluster mass function and astrophysical 
studies of cluster formation. Furthermore this method will probably be the least expensive way to 
measure masses of large numbers of z > 1 clusters detected in future IR-imaging surveys. 
Subject headings: galaxies: photometry 



1. INTRODUCTION 

Galaxy clusters represent the largest collapsed struc- 
tures in t he Universe which form in the peaks of the den- 
sity field (iVoitI 120051 1 Allen et al.ll201lL for reviews see). 
In order to relate the observed cluster population to the 
density peaks seen in numerical simulations the masses 
of the observed clusters have to be inferred from observ- 
ables. An accurate calibration of these mass-observable 
relations is of great importance in astrophysical studies 
of cluster-formation and -evolution as well as in cosmo- 
logical studies using the galaxy cluster ma ss func tion 

The richnes s of a ga laxy cluster (lAbelll 

Postman et al.l 119961: lYee fc L6pez-Cru3 
Muzzin et aLri2007al: iRozo et al.l I2009D is such 
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observable that can easily be measured from optical/IR 
data out to high redshifts. The mass measurements 
that are required for calibration, however, become 
increasingly difficult at higher redshifts. The strong 
assumptions that are necessary to interpret X-ray 
(hydrostatic equilibrium) or velocity dispersion (virial 
equilibrium) measurements make it important to per- 
form independent cross-checks which do not rely on the 
same assumptions. This is particularly important at 
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high-z where clusters are dynamically younger and equi- 
librium assumptions become increasingly questionable. 

One way to measure mass that does not rely on as- 
sumptions like hydrostatic or virial equilibrium is weak 
gra vitational lensing (WL; for a comprehensive review 
see lBartelmann fc Schneideiil200ll) . The commonly used 
measurement of the shearing of galaxies breaks down 
for higher redshifts because the sources cannot be re- 
solved anymore and hence their ellipticities cannot be 
measured. This is particularly true for ground-based 
data and lens redshifts approaching z ^ 1. Those higher 
redshift regions are, however, of special interest for clus- 
ter astrophysics because they approach the cluster for- 
mation epoch where observational guidance is crucial. 
The measu rements of dark energy using clusters (see e.g. 
iVoiti 120051 ) as we ll as measurements of primordial non- 
Gaussianity fe.g. iMortonson et al.ll20lil ) benefit from a 
large redshift baseline including high-z observations. 

Here we present a first measurement of the average 
mass of such high-z cluster lenses b ased on the magnifi- 
cation effect of WL as proposed by Ivan Waerbeke et al.l 
(|2010D . This method does not require one to resolve 
the sources and hence can be applied at higher red- 
shifts. By stacking a number of clusters we can mea- 
sure their average mass. The clusters are selected 
from the Spitzer Adaptatio n of the Red-sequence Clus- 
ter S u rvey (SpARC jH; Wils on et al.ll2009l: IMuzzin et all 
120091 : iDemarco et al.i i2010il : whereas the background 
sources are selected from catalogues of the Canada 
France Hawaii Telescope Lensing Survey (CFHTLenS), 
which are based on imaging data from the CFHT 
Legacy Survey (CFH TLS) . via the Lyman -break tech- 
nique ()Steidel et al.lll996t iGiavaliscol [20021) . The data 
are presented in Sect. [H In Sect. [3] we describe the 

1 |http : //faculty ■ ucr ■ edu/$\slm$gllllanw/SpARCS | 
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magnification method. Results are reported in Sect. 2] 
and discussed in Sect. [5j Throughout we assume a flat 
ACDM cosmology with D,^ = 0.27, VI a = 0.73, and 
Hp = 70km~^s~^Mpc^^ according to iKomatsu et al.l 
()2011|) . All magnitudes are in the AB system. 

2. DATA SET 

The cluster sample used for this analysis is drawn 
from the 9 deg^ XMM-LSS field which is one of six 
fields in the 42 deg^ SpARCS survey. Clusters are se- 
lected as overdensities in a combined color and position 
sp ace using the clu s ter red-sequ ence method developed 
by IGladders fc Ye^ (l2?)00l [2005h . SpARCS uses a z/ - 
3.6/j,m color selection which allows the detection of clus- 
ters up to redshifts as high as z ~ 1.5. Photometric 
redshifts for the clusters are calculated based on the 
color of the cluster red-sequence and have been refined 
using a subsample of 33 clusters with confirmed spec- 
troscopic redshifts. Comparison with the spectroscopic 
sample shows that the photometric redshifts are accurate 
to Az 0.1 up to z ~ 1.4. The richness is estimated from 
3.6/im data. Since the 3.6/im color of galaxies is inde- 
pendent of star formation history and redshift for z <2, 
these 3.6/im richness estimates are equivalent to K-hand 
ones. Further d etails of the S pARC S cl uster catalogs 
can b e found in iMuzzin et al.l (|2009f ) and iWilson et al.l 
(|2nn9t) . and a more detailed discussion of the cluster de- 
tection algorithm can be found in Muzzin et al. (2008). 

In the following we concentrate on secure clusters 
(SpARCS fiux > 9). This ensures that only massive sys- 
tems and no potential false positives are stacked. Select- 
ing clusters in this way leads to more massive systems at 
high redshift than at low redshift, but with the limited 
data set used here we cannot afford more physical cuts 
(e.g. fixed richness). Within the overlap area between 
SpARCS and CFHTLS we find 48 moderately rich and 
securely detected clusters with redshifts between z = 0.2 
and z = 1.6. 

As sources we use Lyman-break galaxies (LBGs) se- 
lected from catalogues of the CFHTLS field Wl which 
overlaps with the SpARCS XMM-LSS field by ap- 
proximately six deg^. The optical mult i-color cata- 
logue s are created from stacked images (see lErben et al.l 
120091 for details on the imaging data reduction with 
the THELI pipeli n e) in a similar way as described in 
iHildebran dt et ahl (|2009a| ) with the notable exception 
that the PSF is brought to the same Gaussian shape over 
the whole pointing and for all five images in the ugriz- 
filters of each field. This allows more accurate color mea- 
surements which results in better photo-z's and a more 
accurate Lyman-break galaxy selection. Details of this 
method will be presented in a forthcoming paper. The 
LEG color selection of u -, q-, and r-dropou t s is ide ntical 
to the one presented in IHildebrandt et ahl ()2009aD with 
the exception that u-dropouts are only required to have 
u — g > 1 instead of it — g > 1.5. This relaxed u — q 
cut is identical to the one used by Ivan der Burg et all 
()2010t ) whose luminosity function measurements we use 
for calibrations 

In the following we use z ^ 3 w-dropouts and z ~ 4 
g-dropouts. The u-dropouts are the cleanest background 

^ This change was introduced to ma ke the samp l e mor e compa- 
rable to other u-dropout studies, like .Steidel et aLl l|1999l ). 



sample in terms of low-z contamination and they are 
the lowest-redshift dropout sample detectable from the 
ground which results in the highest apparent brightness. 
The (/-dropouts are fainter by ^ 0.5mag due to their 
larger distance resulting in larger photometric errors. 
They also show larger low-z contamination. Thus, we 
can only safely cross-correlate them to clusters with in- 
termediate and high redshifts - beyond the redshifts of 
the possible contaminants. In the following, we correct 
for the dilution due to contamination (boosting the signal 
by 10%) whenever g-dropouts are used. The CFHTLS- 
Wide data are too shallow to select r-dropouts in signif- 
icant numbers. 

3. MASS MEASUREMENT 

The signal-to-noise ratio (S/N) of the lensing signal per 
background galaxy is generally lower for magnification- 
based methods than for shear-based ones|3 However, 
since magnitudes - the only requirement for magnifica- 
tion - are easier to measure than shapes, there are always 
more galaxies and in particular higher-redshift galaxies 
available for magnification. Thus, there is a break-even 
redshift for each data set beyond which magnification 
becomes more powerful than shear because ellipticity 
measurements become impossible. This is the case with 
clusters at z > 0.8 where very few galaxies with reli- 
able shape measurements are available at redshift higher 
than the cluster redshift, as it is very difficult to measure 
shapes of galaxies with z > 1 from ground-based data. 

Since the S /N is too low for a single high-z cluster to be 
detected via magnification with the background samples 
described above we rely on stacking the signals of several 
clusters. In this way we can in principle estimate their 
average mass, a method very similar to what is usually 
done in galaxy-galaxy-le nsing. Therefore, the estim ator 
is identical to Eq. 11 in IHildebrandt et al.l ()2009b[ ). i.e. 
we look for correlations in the positions of LBGs and the 
SpARCS clusters. Those positions are correlated due to 
the magnification-bias effect of WL. 

The magnification signal scales linearly with the slope 
of the number counts, a|3 Thus, it is important to select 
a source population with a large |a|. For the interpre- 
tation of the signal, a needs to be measured. However, 
we cannot use the measured a directly because of incom- 
pleteness in our catalogue. Furthermore, the incomplete- 
ness changes over the magnitude bin and so does the a 
itself. In order to account for these effects we model the 
incompleteness as a function of magnitude by comparing 
the LBG number counts in the CFHTLS- Wide to th e 
ones in the CFHTLS-Deep (|Hildebrandt et al.l [2009af ) . 
We then m ultiply the luminosity fun ction of the LBGs 
measured in Ivan der Burg etall poTo.f l with this incom- 
pleteness function. Next we artificially magnify the same 
luminosity functions for different values of the magnifica- 
tion /X and again multiply with the un-magnified incom- 

^ This depends on th e slope of the number counts (see e.g. 
IHildebrandt et al.ll2009bh . It is beneficial to select a source popu- 
lation with a steep slope. 

^ The a we use here should not be confused with the faint-end 
slope of the luminosity function which is also often called a. 

^ Note that the volume probed by each of the dropout samples in 
the four deg^ of the CFHTLS-Deep fields is ~ 0.2Gpc^/i-i so that 
cosmic variance on the number counts is negligible. For compari- 
son: the Millennium Simulation has a volume of 0.125Gpc''h~^ . 
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pleteness function. Then we compare the numbers of 
LBGs in a chosen magnitude bin between the un-lensed 
and the lensed case. A linear fit to the number excess 
(or depletion) as a function of /i yields Ofcff, the effective 
a to be used in the interpretation of the measurement. 

The two problems of varying incompleteness and vary- 
ing a over the magnitude bin should be considered in all 
magnification applications. Just correcting the number 
counts for incompleteness before measuring a and then 
averaging the individual a's of all galaxies (assigned de- 
pending on their magnitudes) in the magnitude bin only 
works if the changes of the incompleteness as well as of a 
over the magnitude range are moderate. The aeff's that 
we use here can easily differ from these naive a estimates 
by 20-30%. 

Besides estimating accurate aeff's it is important to 
cleanly separate the lenses and the sources in redshift 
to prevent physical cross-correlations from biasing the 
WL measurement. Using LBGs as a background sample 
(and M-dropouts in particular) has the great advantage 
that they represent a very clean high-z sample and their 
small amount of low-redshift contaminants are very well 
understood (Steidel et al. 1999, 2003). In this way we 
can make sure that no galaxies with redshifts comparable 
to the cluster redshifts are included in the background 
samples and that physical cross-correlation (in contrast 
to magnification-induced cross-correlations) do not play 
a role here. 

The cross-correlation signal wjd) is then related to 
the m agnification fi through ([Bartelmann fc Schneideil 
I2OOI : 

w{e) = (aeff - 1) s^ii9) , (1) 

with 6fi{e) = pL{e) - 1. The (acff - 1) for the bright 
LBGs selected from the CFHTLS-Wide, which we use 
here, is always positive, so that we expect a positive 
cross-correlation signal. 

Since the clusters span a range of redshifts and masses 
- as suggested by their richness - we cannot just fit a sim- 
ple cluster mass model like a singular-isothermal-sphere 
(SIS) to the stacked signal and estimate the average clus- 
ter mass directly. Rather we fit the WL signal with a one- 
parameter model which is constructed from a weighted 
sum of SIS models for all N clusters: 

N 

Mcomb(a) = fisisjzi; a M200 .richness, (2) 

1=0 

The relative weighting between the clusters in this 
model takes the different redshifts Zi (and hence different 
lensing efficiencies and different critical densities) as well 
as the different masses M2no.richness.z1 a. s estima ted from 
their richness via Eq. 9 of IM 

uzzm etan (I2007bn . into ac- 
count. Hence the fit does not directly yield an average 
mass. Instead the normalization of the mass-richness re- 
lation is fitted under the assumptions that the slope of 
this relation does not change with redshift, the redshifts 
for the clusters are accurate, and that the background 
cosmology is known. 
The magnification for an SIS takes the form: 



with 9e = 47r (^) being the Einstein radius of an 



SIS with velocity dispersion ct„ ; Dg and D^s are the angu- 
lar diameter distances from observer to source and from 
deflector to source, respectively. The velocity dispersion 
is related to the mass, M200, by 



''v = ^M2oo]l -^200pcnt(^)G3 , (4) 

where Pcrit(-z) is the critical density at the cluster red- 
shift, z, and G is the gravitational constant. Here we 
estimate the masses with the mass -richness relation at 
z ^ 0.3 from iMuzzin et al.l ()2007bO which is based on 
dynamical mass estimates. A fitted a of unity would 
mean that the mass-richness relation at t he redshifts of 
our clu ster is identical to the one from IMuzzin et al.l 
(l2007bD . The mass estimated from magnification be- 
comes ^^magnification — ^ -M^j-jchncss • 

In the future, larger cluster samples will allow us to bin 
the clusters much more finely so that an equal weighting 
can be used and no external mass-observable relation is 
needed for proper scaling. 

4. RESULTS 

Fig. [T] shows the cross-correlation signal between rich 
SpARCS clusters in two different redshift bins and the u- 
dropouts with 23 < r < 24.5 ((acff — 1) = 1-54; surface 
density 0.26 arcmin"^) and g-dropouts with 23.5 < i < 
25 ((ofeff - 1) = 1.24; surface density 0.24 arcmin^^)13 
The mean masses of the clusters are estimated from 
their richness. Errors on the data points represent the 
standard deviation between the signals for the differ- 
ent clusters and thus reflect Poissonian noise as well as 
noise due to clustering of the background population and 
noise introd uced by structures alo ng the lines-of-sight to 
the clusters (jHoekstra et al.l[2011[ ). The off-diagonal ele- 
ments of the covariance matrix are estimated by measur- 
ing the cross-correlation signal between 1000 fake clusters 
placed at random positions and the background LBGs. 

We fit Multi-SIS (sec Eq. ^ profiles to the combined 
u- and (/-dropout signals for the two cluster samples in 
the range 0.'8 < 9 < 20'. Table [T] summarizes the results. 
The WL signal is - within the la rge errors - consistent 
with the signal expected from the IMuzzin et al.l ()2007b[ ) 
mass-richness relation, as indicated by the values for a 
being consistent with unity. Despite similar S/N the rela- 
tive error of the mass is larger for the high-z sample than 
for the low-z sample. That is due to the fact that the 
high-z g-dropout signal is not very well-fit by the multi- 
SIS profile (but still consistent within noise). Note that 
we just used the scaling with richness to construct the 
multi-SIS function and not the absolute normalization of 
the mass-richness relation. Separate multi-SIS fits to the 
signals of the u- and g-dropouts of a given cluster sam- 
ple yield consistent results as well. In the following we 
perform two tests for systematic errors that could affect 
our result. 

The SpARCS XMM-LSS field also overlaps with the 
CFHTLS-Deep field Dl. In this one deg^ field much 
deeper imaging is available which allows the selection of 

® When comparing surface densities between shear and magnifi- 
cation to estimate the S/N, the magnification densities should be 
scaled by a factor (2|a(,ff ~ ll^e) ijSchneider et al]|2000l 1. with (Te 
being the intrinsic ellipticity dispersion. 
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Figure 1. Magnification-induced cross-correlation between rich SpARCS cluster and u-dropouts (top) and g-dropouts (bottom) from 
CFHTLenS. The dotted lines correspond to the best-fit multi-SIS profiles and the S/N values are for the individual measurements (see 
Table [T] for the S/N when u- and g-dropouts are combined). Left: All 48 clusters. Middle: Subsample of 32 clusters with z < 0.5. Right: 
Subsample 16 clusters with z > 0.5. 



Table 1 

Results of the multi-SIS fits to the magnification signals. 
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fainter LBGs (the ones also used in iHildebrandt et al] 
l2009ai rb[) which have negative {aes — 1)- We find eight 
of the 48 clusters in that field. Although the statistical 
power is limited due to the small number of clusters and 
the lower absolute value of {aeS — 1) (somewhat com- 
pensated by the higher LBG density), these faint LBGs 
allow us to perform an important test. While the positive 
cross-correlations of the brighter LBGs could in principle 
also be caused by improper redshift separation and hence 
physical cross-correlations, anti-correlations can only be 
caused by WL. Here we use faint u- and g-dropouts 
with r > 25.8 and i > 26.3, respectively, which corre- 
spond to parts of the LBG luminosity function where 
(aoff' - 1) = -0.09 and {aeS - 1) = -0.18. Their surface 
densities are 2.7arcmin~^ and 0.83 arcmin"^. 

As expected for a background sample with {aoB — 1) < 
0, we find under-densities of LBGs around the clusters. 
The detection is less significant due to the limited num- 
ber of clusters and the lower absolute values of (aeS — 1) 
which cannot be compensated by the higher surface den- 
sity of dropouts in the CFHTLS-Deep and their lower 
intrinsic clustering signalQ 

Another useful test is to cross-correlate the clusters to 
a population of objects that should not show any correla- 
tion in their positions because it is closer to the observer 
than the clusters. For that test we select stars and mea- 
sure their cross-correlation function to the clusters. We 

''' The main noise term in the magnification measurement is not 
the Poissonian shot-noise, which originates from the finite sampling 
of the K-field with a limited number of background objects, but 
the noise int roduced by the intri nsic clustering of the background 
galaxies (see Ivan Waerbeke"2010'. for a full treatment of the differ- 
ent error terms) . As we showed in Hildebrandt et al. (2007, 2009^) 
the clustering of the faint LBGs is much weaker than the clustering 
of the bright LBGs. 



do not detect any significant cross-correlation signal for 
stellar samples of different limiting magnitudes. It should 
be noted that the stellar auto-correlation function is dif- 
ferent from zero (with a very low amplitude, but high sig- 
nificance) so that a non-detection of a cross-correlation 
signal between stars and clusters is a meaningful test. 

From these two tests we conclude that our measure- 
ment is robust and free of strong systematic biases. How- 
ever, the current sample is two small to study more subtle 
systematic effects in detail. 

5. DISCUSSION AND OUTLOOK 

We measure the normalization of the high-redshift 
mass-richness relation of galaxy clusters selected from 
SpARCS employing a novel method based on the WL 
magnification effect of these clusters on z > 3 LBGs. 
Even with a small sample of rich, high-z clusters and 
optical data of modest depths we detect a magnifica- 
tion signal at > 3-cr significance which yields a mea- 
surement of this normalization (equivalent to a measure- 
ment of the average cluster mass) with a relative error 
of ~ 50% The measured normalization is in good 
agreement with the low-z mass-richness relation from 
Im uzzin et aTl (|2007b( ) which is based on dynamical mass 
estimates. 

Using much deeper data on a one deg^ sub-field we de- 
tect an anti-correlation between a few high-z SpARCS 
clusters and faint LBGs supporting a WL origin of the 
signal and suggesting that the measurement on the shal- 
lower data is free of physical cross-correlations that could 
mimic a WL signal. 

The overlap area between SpARCS and CFHTLS is 
too small to subdivide the clusters into different rich- 
ness bins and estimate a full mass-richness relation at 
high redshift. However, it is well within reach of current 
facilities to expand this area by an order of magnitude 
through either observing IR-, X-ray-, or SZ-surveys with 
ground-based optical imaging or through targeting op- 
tical surveys such as the CFHTLS with IR- or X-ray 
imaging from space or with SZ telescopes. 



* Note that the relative error of a scales as 
that agree within errors with the model. 
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We would like to stress that this method represents a 
way to obtain mass estimates of large samples of clus- 
ters at very high redshifts that is relatively cheap in 
terms of telescope time. It is the only method that can 
realistically be applied to measure the masses of very 
large numbers of 2; > 1 clusters to be found by future 
high-z cluster surveys. Concentrating on the bright part 
of the LBG luminosity function and probably also on 
the lowest redshift and hence brightest LBGs reachable 
from the ground (w-dropouts) seems the most promising 
route. This certainly requires high-quality u-band imag- 
ing which is available at a few facilities only. 

Based on observations obtained with 
MegaPrime/MegaCam, a joint project of CFHT 
and CEA/DAPNIA, at the CFHT which is operated 
by the NRC of Canada, the CNRS of France, and 
the University of Hawaii and on observations made 
with the Spitzer Space Telescope, which is operated by 
the Jet Propulsion Laboratory, California Institute of 
Technology under a contract with NASA. This work is 
based in part on data products produced at TERAPIX 
and CADC. 
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III" and by the DFG through project ER 327/3-1 and 
the TR 33. H. Hoekstra acknowledges support from 
the NWO and a Marie Curie IRG. LVW is supported 
by NSERC and ClfAR. GW acknowledges support from 
NSF grant AST-0909198. 
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